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Abstract 

We consider pion-nucleon scattering in heavy-baryon chiral perturbation theory to third 
order. All electromagnetic corrections appearing to this order are included. We thus have 
a consistent description of strong and electromagnetic effects, which allows us to isolate 
the strong part of the interaction in an unambiguous way. We give pion-nucleon phase 
shifts up to pion laboratory momenta of 100 MeV and find sizeable differences for the 
S-waves of the elastic channels compared to previous phase shift analyses. The precise 
description of the scattering process also allows us to address the question of isospin 
violation in the strong interaction. For the usually employed triangle relation we find 
an isospin breaking effect of —0.7% in the S-wave, whereas the P-waves show an effect 
of —1.5% and —4 to —2.5%, respectively, for pion laboratory momenta between 25 and 
100 MeV. 
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1 Introduction 



The study of the strong interactions at low energies has stimulated vast efforts in the domain 
of pion-nucleon physics. It is mainly the description of simple processes like pion-nucleon 
scattering which can help us to extend our knowledge about the fundamental principles of the 
strong interactions. For every field theory, the consideration of symmetry principles is very 
important, and it is thus of primordial interest to know to what extent these symmetries are 
realized in nature. 

For a long time, isospin symmetry has been thought of as an exact symmetry of the strong 
interactions. But the nucleon isospin doublet, consisting of proton and neutron, is not degen- 
erate in mass; whereas electromagnetic interactions alone would make the proton heavier than 
the neutron, the mass difference between the up- and down-quarks reverses the situation: the 
neutron turns out to be heavier than the proton by 1.3 MeV. This simple consideration triggers 
questions about the size of isospin violation in other strong processes. 

Pion-nucleon scattering is the ideal playground for studying the question of isospin breaking. 
Indeed, thanks to the construction of meson factories and the formidable effort of many experi- 
mental groups, we now have a large amount of very precise low-energy data. Weinberg pointed 
out on purely theoretical grounds that isospin breaking effects in neutral-pion scattering must 
be dramatically enhanced due to the smallness of the isoscalar amplitudes |J. The advent of 
high-quality data opened the way to the study of isospin symmetry violation in processes which 
are experimentally accessible: these are the elastic ir + p and n~p as well as the single charge 
exchange (SCX) reactions. The difficulty now lies in the extraction of the strong interaction 
amplitude from the scattering data. Electromagnetic corrections have to be applied in order to 
gain information about the hadronic amplitudes. Such corrections have been computed e.g. by 
the Nordita group and more recently by Gashi et al. [[3] . The problem is that these analyses 
cannot account for the electromagnetic and the strong interaction in a consistent way. In view 
of studying small quantities like isospin symmetry violation, it is absolutely necessary to avoid 
systematic errors resulting from a possible mismatch in the description of both forces. We will 
readdress this question, i.e. we will extract the strong amplitude from scattering data in the 
framework of chiral perturbation theory (%PT). ;\PT allows to calculate the electromagnetic 
and the strong interactions consistently in terms of unknown parameters, the so-called low- 
energy constants (LECs). Having fitted the full amplitude to the available cross section data, 
one is then able to switch off electromagnetism and to predict the strong piece in a unique 
fashion. 

The knowledge of the hadronic amplitude is necessary for analyzing the violation of isospin sym- 
metry of the strong interaction. Such analyses have been performed in the past ||, [5] , indicating 
breaking at low energies as large as 6%. Both works, however, are based on electromagnetic 
corrections which might not be compatible with the description of the strong interactions. In 
order to better understand the origin of this large isospin breaking effect, we readdress this 
question within our framework. The advantage of such a diagrammatic calculation lies in the 
fact that we can clearly localize which mechanisms are responsible for isospin violation, thus 
making a detailed study of the subject feasible. 
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2 The amplitude 



This work follows a series of previous papers: in ||, pion-nucleon scattering was studied in 
the limit of isospin symmetry. The unknown counterterms were determined by a fit to phase 



shifts resulting from three different analyses ]7|, g, ||. In [[Li], [Tl| , we included strong isospin 
breaking terms in the amplitude and took the full mass difference of the physical particles into 
account. The isospin symmetric amplitude had to be generalized in order to describe scattering 
of particles with different masses. Not only do the particles in the initial and final states have 
different masses, but pion mass differences also have to be considered in the loops, and the 
nucleon mass difference plays a role in the intermediate states. 

The aim of this work is to describe QCD and QED with unequal up- and down-quark masses 
(it^u 7^ f^d) and non-vanishing electric charge (e 2 7^ 0). We thus have to add loops with virtual 
photons as well as electromagnetic counterterms. This is the yet missing part in the complete 
analysis of pion-nucleon scattering. Photons enter the amplitude via four mechanisms, one 
example of each is shown in fig. [I]. 

• First, there is the Coulomb potential, which is the exchange of soft photons between the 
charged pions and the proton. To the accuracy we are working we only have to include the 
one- and two-photon exchange, since multi-photon exchange diagrams are suppressed by 
higher powers of e 2 . It is important to note that our simultaneous counting of momenta 
q and electromagnetic couplings e as small quantities necessarily breaks down at very 
low energies. To illustrate this, let us compare the Weinberg-Tomozawa term (which 
in our counting is of first order) and the same process accompanied by the exchange of 
a photon between the initial or the final particles. At low energies, the loop diagram 
(which is of third order) is "suppressed" by a factor e 2 M ff /(32|g^.|). Thus, for pions of 
momentum less than 0.5 MeV, the loop diagram is not suppressed compared to the tree 
diagram, which clearly contradicts our counting scheme. At very low energies, one would 
thus have to sum up diagrams with an arbitrary number of photons exchanged between 
the charged-particle legs. This is commonly done in nucleon-nucleon scattering, where 
the heavy nucleons can be considered to be non-relativistic at such energies. Fortunately 
we do not have to worry about these ladder exchange diagrams since they only become 
important at energies much lower than the experimentally accessible energy range. In our 
analysis, we will thus stick to the traditional chiral counting, which consists in considering 
e a quantity of chiral order q. 

• Second, there are the "usual" loop diagrams with virtual photons, i.e. those which include 
at least one strong-interaction vertex. Many of these only renormalize the coupling 
constants. 

• Third, we also have to consider Bremsstrahlung processes where a soft photon is emitted 
from an external particle, but is not detected since its energy is lower than the detector 
resolution AE. The inclusion of this process makes the amplitude infrared finite. We have 
not performed this calculation analytically, since it is sufficient for the actual computation 
to combine every log(m 7 /M 7r ) from virtual-photon diagrams with \og{AE/m 1 ) due to the 
integrated emission of soft photons with energy smaller than the detector resolution. The 
divergence appearing because of the vanishing photon mass m 7 thus cancels out. For the 
evaluation, we choose AE = 10 MeVP 3 ]. 



# 3 This choice is rather arbitrary; we have checked, however, that our results do not depend on the particular 
choice of AE 
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• Fourth, hard photons have been integrated out of the theory. Their effect is parameter- 
ized by a string of terms in the Lagrangian, each accompanied by an unknown coupling 
constant. These constants, similarly to the hadronic counterterms, have to be fitted to 
experimental data. 




a) b) c) d) 

Figure 1: Example diagrams for a) the Coulomb potential, b) a "usual" photon loop, 
c) Bremsstrahlung, and d) an electromagnetic counterterm. The solid, dashed, and wiggly 
lines denote nucleons, pions, and photons, in order. 



The inclusion of photons and electromagnetic counterterms in the amplitude gives us a de- 
scription of pion-nucleon scattering based on full QCD and QED. Lepton loops, e.g. vacuum 
polarization, have not been taken into account explicitly, but their effect is hidden in the elec- 
tromagnetic counterterms. One could equally well compute lepton loops by using a Lagrangian 
containing photons and leptons as active degrees of freedom fll|, but for our purposes the 



parameterization of these effects by counterterms is sufficient. Let us add a last remark about 
the relativistic form of our amplitude. The pions and photons are treated fully relativistically; 
on the other hand, since we restrict ourselves to the energy region just above threshold, we 
can consider the nucleons to be very heavy and work in heavy-baryon chiral perturbation the- 



ory [13|, |14 |. This is a particular non-relativistic expansion of the amplitude, which amounts 
to organizing observables in a series of terms with increasing power of the inverse (heavy-) 
nucleon mass. 

Our calculations are based on the Lagrangian constructed in [15| (see also |16|)- The relevant 
terms can be found in appendix 0. 



3 Fitting procedure 

After having constructed the full amplitude, we can proceed to fixing the unknown parameters. 
For this we use cross section data for elastic n + p and 7t~p scattering as well as for the charge 
exchange 7i~p — > ir°n reaction. Guided by the results of [fj, we decide to use experimental 
information up to pion laboratory momenta of 100 MeV (i.e. T n < 32.2 MeV) as input. We 
perform a combined fit of our isospin breaking amplitude to all three channels which have been 



studied experimentally. In the elastic channels, we use the data by Frank et al. [17|, Brack 



et al. [ITS]], and by Joram et al. Low-energetic charge exchange data have been taken by 
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Isenhower et al. [^0|, Salomon et al. 
the old data by Bertin et al 
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13], and by Frlez et al. |22 
and by Duclos et al. 
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Note that we do not include 
in our data base, since these are 



nowadays generally believed to be erroneous. 

In the full amplitude, there appear 14 combinations of parameters which have to be fixed. Such 
a large number of unknowns makes a determination from pion-nucleon scattering data alone 
(almost) impossible. But, fortunately, some of these counterterms appear in other processes, 
which have been studied in great detail. The counterterms c 6 /7, <i 6 / 7 , related to the coupling of 
a photon to the nucleon, can be fitted to the magnetic moments and the charge radii of proton 
and neutron (see e.g. [25|). The separation of the nucleon mass difference into a strong and an 
electromagnetic piece [^B], [^J fixes C5 and fi- We are now left with 11 parameters which have 
to be determined. 



We define the function which has to be minimized in analogy to the one introduced in [28 



X 



with 



Xi 



j 



(.Vij Vij 



stat N 2 
ij I 



the partial contribution of the j'th data set comprising rij 
given by 



j 



i.e. 

(1) 

(2) 

entries, and the weights Wj being 

(3) 



Here N is the number of data sets we fit to, and the uncertainties are split into a statistical piece 
<rf* at and a systematic contribution <r^ ys 0. The experimental data are denoted by ?/^ xp , and we 



choose the LECs in such a way that our predicted cross sections yfj minimize the x 2 ~function. 

We use the standard MINUIT minimization routines of the CERN library. Table |l] gives our 
best-fit values at the scale A = M n o. Note that the quoted errors correspond to the MINUIT 
errors only and are certainly underestimated. The hadronic low-energy constants are of natural 
size, i.e. of order 1; as for the electromagnetic counterterms, some come out unnaturally large, 
but in this case they are accompanied by big uncertainties. It is obvious that one has to be 
careful about correlated parameters in a fit involving such a large number of unknowns; we 
indeed have large correlations among the strong-interaction counterterms. This does not come 
as a surprise since we restrict our analysis to a small energy region above threshold. There 
the pion energy, for example, is hardly distinguishable from the pion mass which causes the 
parameters accompanying these two structures to be strongly correlated. For these reasons 
we refrain from making predictions for quantities like the sigma term or the pion-nucleon 
coupling constant g^NN which are determined by LEC combinations strongly correlated to 
other parameters. We have checked, however, that there are very small correlations between 
electromagnetic and strong counterterms^ 5 !. This guarantees that we can reliably extract the 
strong phases from the full amplitude. The smallness of the x 2 /dof value shows that we fit to 
data in a region where the third order amplitude is still precise enough. More important, it 
also demonstrates the consistency of the low-energy data base. 

^ 4 (T^ ys is determined by the normalization uncertainty of the specific experiment. In analogy to ]28| ], we 
attribute a systematic uncertainty of 5% to data sets for which this quantity is not given explicitly. 
* 5 ci needs special treatment in this respect. We will come back to this point in section H. 
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Table 1: Values of the LECs resulting from a fit of the full pion-nucleon amplitude to low- 
energy cross section data. Parameter values marked by a * are input quantities. 



4 Extraction of the strong amplitude 

After having determined all low-energy constants appearing in the pion-nucleon amplitude, we 
can now proceed to extracting the strong interaction piece. As pointed out in the introduction, 
it is this piece which is of particular interest for the determination of isospin breaking. In the 
previous paragraph we have used the full amplitude in order to describe pion-nucleon cross 
sections. The strong amplitude is defined as the QCD scattering amplitude with m u ^ rrid, but 
e 2 = 0. We will thus not work in the isospin symmetry limit, but we will leave room for strong 
isospin breaking. 

In order to determine the strong amplitude, we have to know what the hadronic masses of the 
scattering particles are. We can write the charged and neutral pion masses as 

M 2 + =Mt+^e 2 C+Afe* + Bfe 2 {m u -m d )+C ( t\m u -m d ) 2 , (4) 

M 2 = M 2 +4 4) e 4 + Bfe 2 [m u - m d ) + C<) 4) (m« - m d ) 2 . (5) 

Here Mo is the bare pion mass (more precisely, the leading term in the quark mass expansion 
of the pion mass in the two-flavor case with m s fixed), and the «4+/ , £>+/ and C+j Q are 
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combinations of low-energy constants from the fourth order mesonic Lagrangian with virtual 
photons. Note that the pion mass difference to leading order is of purely electromagnetic origin. 
We obtain the hadronic pion masses by switching off the electromagnetic interaction and using 



\MZo-M n o\ 
ML 



< 1 MeV , 

= MK + 0.1 MeV 



(6) 
(7) 



These numbers are based on the work by Knecht and Urech |29| using naturalness of the 
low-energy constants. 



As for the nucleon masses, we have 



m p /n = m- AMloCi =F 2B (m u - m d )c 5 - e 2 F 2 (7i + + fs) 



(8) 



in terms of the bare nucleon mass m (i.e. the nucleon mass in the SU(2) chiral limit) and the 
counterterms given in appendix |A|. The hadronic masses will then be given by 



m - AMloCi - 2B (m u - m d )c 5 
m p + e 2 F 2 (f 1 + f 2 + f 3 ) , 
m - AMloCi + 2B (m u - m d )c 5 
m n + e 2 Fl{h + h) . 



(9) 
(10) 



In the following we will take these mass differences into account, as well for the external particles 
as for the nucleons and pions in the intermediate statesp 6 . 



It becomes clear from eqs. (§) and (10) that we will not be able to determine the hadronic 
nucleon masses precisely, since we do not know the value of the LEC combination (fi + / 3 ). 
An additional problem is due to the fact that, after including electromagnetic counterterms 
into the amplitude, we can no more determine C\ and fi separately, as only the combination 
ci = c\ + e 2 F 2 fi/(2M 2 () ) appears. Since we do not know the precise value of these LECs, all 
we can do is to assume natural sizes for f\ and f 3 ; we will thus vary these parameters over 
the interval /1/3 = (Oil) GeV -1 [II|. Our results will, however, turn out to be insensitive 



to the precise values of these parameters. Also the lack of knowledge of the absolute value of 
the hadronic pion masses (see eq. @) does not affect our results; only the mass differences are 
relevant. 

Fig. H shows our prediction for the strong phases as compared to the phase shift analyses 
performed by Matsinos et al. 0, |28| which do not contain isospin violating contributions. The 
range of uncertainty in our prediction, resulting from the correlated error bars of the low- 
energy constants, is given by the dashed curves. In the P3-waves, the agreement between both 
analyses is perfect. In the Pi-waves, our solution is closer to the recent analysis EM00. The 
most dramatic discrepancy lies in the S-waves of the elastic channels; this difference is not due 
to the inclusion of the isospin violating effects, but is already present in the isospin symmetric 
part of the amplitude. The kink in the S-wave phase shift for the charge exchange channel 
comes from the fact that the ingoing pion needs to have a lab momentum of ~ 25 MeV in order 

# 6 This has not been done by the previous phase shift analyses which only work in the isospin symmetric limit. 
This is a severe drawback of these analyses, since the inclusion of such effects leads to important changes, as 
will be shown in the following. 
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Figure 2: Strong pion-nucleon phase shifts as a function of the pion laboratory momentum 
q n for the three measured channels. The solid line corresponds to our solution, the dashed one 



to the one-sigma uncertainty range. Also shown are the EM98 || (stars) and the EMOO ]28 
(open squares) phases. 
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Figure 3: Strong pion-nucleon phase shifts as a function of the pion laboratory momentum 
q n for the three measured channels. The solid line corresponds to our solution, the dashed one 
to the one-sigma uncertainty range. Also shown are the KA85 |7]] (full dots) and the SP98 
(open diamonds) phases. 
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to produce the heavier final state. (Remember that only the hadronic mass difference of the 
involved particles is of relevance.) In fig. |3] we show our solution for the phase shifts compared 
to the KA85 |7| and the SP98 |§ analyses. The agreement in the P3-waves is still very good, 
whereas there are some discrepancies in the Pi phases. Our ir + p S-wave phases agree better 
with the KA85 and SP98 ones than with the EM98(00) solution. But there remains a large 
difference in the S-wave 7t~p phase shift. 

In order to explain where the difference in the S-waves comes from, we further analyze the 
contributions to the electromagnetic corrections (see fig. |5]). The strong phase shifts shown in 
figs. |2], [3] are obtained when fitting the full pion-nucleon amplitude with virtual photons to the 
experimental data, and then switching off the electromagnetic coupling. This is also shown by 
the solid line in fig. ^. The dashed line, on the other hand, is obtained if we do not include 
all virtual-photon diagrams, but only the one-photon exchange diagram with dressed vertices 
and the two-photon exchange process (see diagrams a)-c) in fig. ^). The way to treat the 
one-photon and two-photon exchange between the charged pion and the proton is standard. 
The inclusion of corrections to the N'-y and the 7T7 vertex does not lead to big changes in the 
phases. The remaining virtual-loop diagrams (one example of which is drawn in d) in fig. ||) 
have not been included when drawing the dashed line. We think that such diagrams can only 
be accounted for correctly when the strong interaction and electromagnetism are treated in a 
consistent way. These effects are very large and cause this unusual bending of the n~p S-wave 
phase shift. Previous phase shift analyses have obviously underestimated these contributions. 




a) b) c) d) 



Figure 4: a) One-photon exchange diagram, b) two-photon exchange diagram, c) vertex 
corrections to the one-photon exchange, d) typical graph involving non-linear pion-nucleon- 
photon couplings. The solid, dashed, and wiggly lines denote nucleons, pions, and photons, in 
order 



5 Isospin violation of the strong interactions 

The aim of this section is to determine the size of isospin violation in low-energy pion-nucleon 
scattering. Isospin is an approximate symmetry of the strong interactions, for electromagnetism 
it is largely broken. Consequently, it does not make sense to study isospin violation of the full 
pion-nucleon amplitude. The interesting part lies in the strong amplitude, which we have 



9 



s 





50 100 

q„ [MeV] 




-0.2 



50 100 

q„ [MeV] 



Figure 5: Strong pion-nucleon phase shifts as a function of the pion laboratory momentum 
q n . The full line is deduced from an analysis based on the complete pion-nucleon amplitude. 
To obtain the dashed line, we have fitted the strong amplitude together with the one-photon 
exchange diagram with dressed vertices and the two-photon exchange process to the data, and 
then switched off the electromagnetic interaction (see text). 
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isolated in section ^. As in section |], we will switch off every electromagnetic interaction; we 
will thus be left with the isospin violating two-pion-nucleon vertex ~ C5, and the hadronic mass 
differences of the nucleons and the pions. 

In order to quantify isospin breaking, we will employ the usual triangle ratio [|30], [|, [5], I I 



2 p— >tt+ p fir p^ir p V^2 f n p^n°n (11) 

fn+p^n+p fir~p^-n~p V^2 f- n —p~t 7r n 

Here / stands for the strong part of the pion-nucleon amplitude. In the following, we will 
always form the ratio R with the real part of the scattering amplitude^ 7 ]. In case that isospin 
is a good symmetry R vanishes over the whole energy range. 



It has been argued in |TT[ that it is advantageous not to consider the usual P3- and Pi-wave 
projections, but to focus instead on the P-wave parts proportional to the spin-flip and the 
spin-non-flip amplitudes h and g. These are linked to the usual T-matrix elements via 



T 

1 -, 



nN 



I Ei + mi E 2 + m 2 
2m\ V 2m2 



g + ia ■ (q 2 x q\)h . (12) 



E1/2 is the energy of the in-/outgoing nucleon with mass mi/2, the in-/outgoing pions have 
four-momenta qi/ 2 . The usual S- and P-waves are then given by 

S(s) = lQ7T ^ S(s) , (13) 

p s(s) = ^— j= \qi\\q2\\Q{s) + H(s)) , (14) 

y/Ex + m 1 x/E 2 + m 2 ( y 

Pi(s) = iQ^fs m\m\\G{s)-2H{s)) , (15) 



with 



r+l 

S(s) = J ^ dzg(s,z) , (16) 
1 [+ 1 

G{s) = / dzzg(s,z) , (17) 

H{s) = J*dz^±h(s,z) . (18) 

The integration in the last equations is performed over the cosine of the scattering angle z. In 
the following we will consider isospin violation in the S- as well as in the Q— and ^-projections 
ofR. 

There is one additional problem: as pointed out before, the threshold of the charge exchange 
reaction lies at q^ ~ 25 MeV. Below this value, the SCX channel is closed and the triangle 



# 7 In chiral perturbation theory, unitarity is only perturbatively fulfilled. Since imaginary parts are due to 
loop diagrams which only come in at third order, they are known with less precision than the real parts of the 
amplitudes. 
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relation does not make sense. Nevertheless, the different threshold energies of the elastic and 
the inelastic channels are a consequence of strong isospin breaking (remember that only the 
hadronic mass difference of the particles has been included). Therefore, we want to take this 
effect into account; in fig. ^ we show the momentum dependence of R (in %) for momenta 
q n above 25 MeV. In this energy domain all three channels are open. Again, the dashed line 
indicates the one-sigma uncertainty range of our prediction. For the iS-projection, isospin 
breaking lies around —0.7%. Also the ^-projection varies very little with energy, breaking in 
this channel is somewhat larger, ~ —1.5%. The biggest effect is found to be in Ti where isospin 
violation lies between —4 and —2.5% for q n below 100 MeV. 

Instead we could also argue that we do not want to take threshold effects into account, and 
that we want to separate "static" (i.e. related to kinematic effects) from "dynamical" (i.e. due 
to isospin violating pion-nucleon couplings) isospin breaking pH] . In fig. [7], we consider the 
case where the hadronic mass of the proton and the neutron, as well as the hadronic mass of 



the pions, are taken to be equal 



+ m n )/2,ACo = Ml 



There is 



thus no static isospin breaking, but we only show effects due to the two-pion-nucleon vertex 
proportional to C5. This coupling only shows up in the S-wave of the charge exchange reaction; 
it causes isospin to be broken by —0.75% in S. There is no dynamical isospin breaking in the 
P-waves (to the order considered here). 
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Figure 6: R[%] as a function of the pion laboratory momentum q n . The dashed lines indicate 
the one-sigma range of our prediction. 



6 Summary and Outlook 

For a reliable determination of the size of isospin breaking in the strong interactions, it is 
primordial to describe electromagnetic and strong effects consistently. In the present work this 
is achieved by an analysis of pion-nucleon scattering in the framework of chiral perturbation 
theory. ;\PT does not leave any doubt about the correct definition of the hadronic masses 
of pions and nucleons, and allows to extract the strong part of the scattering amplitude in a 
unique way. After determining the unknown low-energy constants by a fit to experimental 
data, we switch off all electromagnetic interactions and describe QCD with unequal up- and 
down-quark masses and e 2 = 0. The so-determined strong phase shifts agree with those of 
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Figure 7: R[%] as a function of the pion laboratory momentum q n in the scenario where 
and M^ = M^+, as described in the text. The dashed lines indicate the one-sigma 
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range of our prediction. 



Matsinos et al. || ^H] in the P-waves, but we find a sizeably different behavior in the S-waves. 
We trace this difference back to the inclusion or omission of non-linear photon-pion-nucleon 
couplings. 

We address the question of isospin violation by studying the usual triangle relation involving 
elastic ir ± p scattering and the charge exchange reaction. An important advantage of the 
calculation lies in the fact that we can easily separate dynamical from static isospin breaking. 
Dynamical isospin breaking only occurs in the S-wave and is very small, ~ 0.75%. Static effects 
do not increase the size of isospin violation in the S-wave significantly; by no means can we 
account for the reported 6 - 7 % isospin breaking |4], ||. 

We have found large error bars on our parameter values. In order to improve this situation, we 
would like to fit to more experimental data. However, a third order %PT calculation allows to 
describe scattering data for pion laboratory momenta not much higher than 100 MeV, a region 
where the data situation is not yet as good as one would hope. A fourth order calculation would 
certainly allow to fit to data higher in energy, but, on the other hand, would also introduce 
many more unknown coupling constants. Since isospin breaking effects are expected to be most 
prominent in the energy region we consider in this work, we do not judge it very promising to 
extend the analysis to full one-loop (fourth) order. Additional data for pion-nucleon scattering 
at very low energies would be very helpful in this respect. Also a combined fit to several 
reactions involving nucleons, pions, and photons, e.g. pion electro- and photoproduction, as 
well as ttN — > ttttN, would help in pinning down the fundamental low-energy constants more 
precisely. 
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A Lagrangian 

In this appendix we give the relevant counterterms, appearing in the pion-pion and pion- 
nucleon Lagrangians. We only display the finite parts. The complete set of terms can be found 



m 



Let us start with the Lagrangian involving pions, nucleons and virtual photons. We only give 
the form of the counterterms, the expression for the terms with fixed coefficients is canonical. 
At second and third order, these read 



C 



(2) 

7*7rA r 



N 



N 



C-2 



C 3 



C 4| 



cx(x+) + u) 1 ) + ^V> + %[S", S"][ii M) u vl 



+hF 2 Q (Q 2 + - Q 2 _) + f 2 F 2 (Q + }Q + + f 3 F 2 (Q + } 2 + ... 



N 



(A.l) 



^-(i(v ■ uujD* + h.c.) + ^-{e^ af3 v a Sp[u^u v }v ■ D + h.c.) 



4m 4m 
^L.(ie>™f>v a S B F+v ■ D + h.c.) - A 



5m z 



ze 



^^F+D. + h.c.) 



-^^(^^{ifj >« ■ D + h.c.) - S^L^ Sa (F+)D fi + h.c 
+z'd 1 [u jtt , [u • D, + z (d 2 - [u^, [D^, u • u]] + ic/ 3 [w ■ u,[v ■ D,v ■ u]] 

+d 5 [X-,v ■ u] + (d 6 - ^) [D", F+W + (d 7 - ^±^) [D», (F+)y 

+d 14 e^v a Sp([v ■ D, + d 15 e^v a S^ Ufl [D u , v u]) 

+g 6 F 2 S»(Q^v ■ u)(Q+) + g 7 F 2 (Q + v ■ u)Q_ 



+g 8 F 2 (Q^vu)Q^ 



N . 



(A.2) 



This heavy-baryon form of the Lagrangian is given in terms of the nucleon velocity and the 
spin operator = \^o [lv v 1 ' . e is the totally antisymmetric tensor in four dimensions, and 
we choose the convention e 0123 = — 1. F is the unrenormalized pion decay constant, N is the 
nucleon spinor, and h.c. stands for the hermitian conjugate of a particular term. Moreover, (A) 
denotes the flavor trace of A, and A = A — (A)/2 its traceless part. 

Let us restrict ourselves to the case in which there are no external sources but where virtual 
photons are included. These are parameterized by the electromagnetic field A^ and the field 
strength tensor F^ v = d^A u — dyA^. The elementary building blocks of the Lagrangian are 
given in terms of the field 



u 



U 



7T Z T ■ 7T 



1/2 



(A.3) 
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and the nucleon charge matrix Q by 

= i{u~ l (d fl + iQA ll )u-u(d ll + iQA ll )u' 1 } , 
X± = u ] xu ] ± ux ] u , 

Q± = -{u^Qu ±uQu' 1 } , 



= -2Q±F^. (A.4) 

X is related to the quark mass matrix, x = 2 _B diag(m n , m rf ), where _B — — (0|^|0)/Fq + . . .. 
We assume Bo Fq, i.e. we work in standard xPT. The covariant derivative reads 

[D„ A] = d,A+[T„A] , 

r M = ^{u\d^ + iQA^u + u(d^ + tQA^} . (A.5) 

In the pion-photon sector, we only need few terms which enter the amplitude via renormaliza- 
tion. In the Feynman gauge, we obtain 

= -\F, u F^ + ^(uX + X + ) + CF^Q 2 + -Q 2 _) , (A.6) 



4*lr = + + (A.7) 



Here, = — |{w^[r M , Q]u ± Q]^} is needed for the renormalization of F w . 
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